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Abstract 

Using data collected by the high energy photoproduction experiment FOCUS 
at Fermilab we study the doubly and singly Cabibbo suppressed decays D + and 
D+ -» K+-K+TT-. Our measurements of T(D + -» K + Tr + Tr-)/T(D + -» K--k + tt + ) = 
0.0065 ± 0.0008 ± 0.0004 and T(D+ -> K+tt+tt')/ F(D+ -» K + K-tt+) = 0.127 ± 
0.007 ± 0.014 are based on samples of 189 ± 24 D + and 567 ± 31 D+ reconstructed 
events, respectively. We also present Dalitz plot analyses of the two decay channels; 
the amplitude analysis of the Df — > K + it + it~ mode is performed for the first time. 



1 Introduction 

A thorough understanding of D hadronic decays requires both branching ratio 
measurements of all the hadronic modes and Dalitz plot analyses to investigate 
the decay dynamics and probe the role of final state interactions. In this 
letter we present branching ratio measurements and Dalitz plot analyses of 
the doubly and singly Cabibbo suppressed decays D + — > K + -k + -k~ and D+ — > 

A" + 7T + 7T~ . 

The naive expectation for the ratio of Doubly Cabibbo Suppressed (DCS) 
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to Cabibbo Favored (CF) branching fractions is tan 4 #c ~ 0.25%. However 
the CF D + — > K~7t + tt + rate could be suppressed by destructive interfer- 
ence between spectator amplitudes with indistinguishable quarks in the final 
state; this is the argument generally proposed to explain the lifetime dif- 
ference between D + and D°. For the hadronic DCS D + decay all the final 
state quarks are different and no destructive interference is present. In this 
simple picture we would expect, neglecting effects of final state interactions, 
t(D+)/t(D°) = T(D CF )/T(D+ F ) = (l/tan 4 ^ c )xr( J D+ cs )/r( J D+ F ). The com- 
parison of the precise FOCUS lifetime ratio [1] t(D + )/t(D°) = 2.538 ± 0.023 
with the branching ratio T(D^ cs )/T(Dq F ) measurement reported here will 
test this interpretation. 

The final D + and samples are selected with cuts to reduce reflections from 
the more copious Cabibbo favored modes and to optimize the signal to noise 
ratios, which are crucial for reliable decay amplitude analyses. Dalitz plot 
analyses have indeed emerged as a unique tool to fully exploit the available 
charm statistics allowing, besides the simple branching ratio evaluation, to 
investigate the underlying decay dynamics. Our understanding of the charm 
decay dynamics has already considerably improved in the last years, but is 
still limited to a few decay channels. The Dalitz plot analysis of the Singly 
Cabibbo Suppressed decay (SCS) D+ — > K + it + it~~ is performed for the first 
time. 



2 Signal selection 

The data for this analysis were collected during the 1996-1997 run of the 
photoproduction experiment FOCUS at Fermilab. The detector, designed and 
used to study the interaction of high-energy photons on a segmented BeO 
target, is a large aperture, fixed-target magnetic spectrometer with excellent 
Cerenkov particle identification and vertexing capabilities. Most of the FOCUS 
experiment and analysis techniques have been described previously [2-4]. 

The suppressed nature of the channels under study requires severe cuts to 
eliminate reflections, both in the signal and sideband regions, from more co- 
pious and favored decays, when one or two charged particles are Cerenkov 
misidentified and/or one neutral particle is missing. Tight cut choices are also 
required to improve the signal-to-noise ratio to perform a reliable Dalitz plot 
analysis. 

The D + and candidates are obtained using a candidate driven vertex 
algorithm. A decay vertex is formed from three reconstructed charged tracks. 
The momentum of the D candidate is used to intersect other reconstructed 
tracks to form a production vertex. The confidence levels (C.L.) of each vertex 
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is required to exceed 1%. From the vertexing algorithm, the variable £, which 
is the separation of the primary and secondary vertices, and its associated 
error ag are calculated. We reduce backgrounds by requiring Ijot > 14 and 10 
for the D + and Df, respectively. The two vertices are also required to satisfy 
isolation conditions. The primary vertex isolation cut requires that a track 
assigned to the decay vertex has a C.L. less than 1% to be included in the 
primary vertex. The secondary vertex isolation cut requires that all remaining 
tracks not assigned to the primary and secondary vertex have a C.L. smaller 
than 0.1% to form a vertex with the D candidate daughters. The decay vertex 
is required to be 3 a outside of the target material to reduce the background 
due to hadronic re-interactions in the target. 

The Cerenkov particle identification is based on likelihood ratios between 
the various stable particle hypotheses [3]. The product of all firing probabil- 
ities for all cells within the three Cerenkov cones produces a x 2 -hke variable 
W / i=-21n(likelihood) where i ranges over the electron, pion, kaon and proton 
hypothesis. For the D + selection, we require AWk = W w — Wk > 4 and 
AW W = Wk — W w > 3.5 for both pions in the final state. For the Df candi- 
dates we require AWk > 4 for the kaon, AW n > 2 for opposite-sign pion and 
AW W > 1 for same-sign pion. These selections minimize the contamination 
from charm background. 

A detailed Monte Carlo study is performed to evaluate the possible contami- 
nations and structures induced by reflections both in the signal and sideband 
regions. 1 The considered sources of background are the D + and three- 
body hadronic decays, where one or two particles are misidentified, and four- 
body hadronic and semileptonic decays, where charged particles are misiden- 
tified and neutrals are missing, namely: D + — > K~it + it + , D+ — > K~K + n + , 
D+ -> K-K + n + , D+ -> 7T+7r-7r+, D+ -> K'tt+ii+u, D+ -> K'K+n+v and 
D + — > K~it + it + -k . For each of them we evaluate the number and the dis- 
tribution shape of events which survive the selection cuts and penetrate into 
the D + and — > K + 7t + tt~ signal and sideband regions. From this study, 
we find that the major contributions for D + come from D + — > K~h + 'k + , 
D+ -> K~K+ir+ and D+ -> X-tt+tt+tt , and for D+ from D + -> K-it+it + , 
D + — > 7r + 7r~7r + and D + — > A' _ 7r + 7r + 7r°. In Fig. 1 the selected K + -n + -n~ com- 
binations for D + and are shown, along with the distribution of the Monte 
Carlo reflected events. We choose to start our fit from the 1.75 GeV/c 2 energy 
threshold since the region below is dominated by the partial reconstruction 
of multi-body charm channels where neutrals are missing. The signal yields 
consist of 189 ± 24 events for D + and 567 ± 31 events for Df. The fits are 
performed with two Gaussian functions for the signals and a second order 



1 The sidebands for D + cover the — 5a to — 3a and the 3a to 5a regions from the 
D + peak. For the the left sideband covers the —5a to —3c region and the right 
sideband covers the 4a to 6a region, both from the Df peak. 
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invariant mass K + n ji + GeV/c 2 invariant mass K + k ji + GeV/o 2 

(a) (b) 

Figure 1. Invariant mass distributions K + 7r + ir~ with final selection criteria for (a) 
D + and (b) Df, along with distributions of the residual reflected events expected 
from Monte Carlo (dotted curves). In the picture the sideband regions (vertical 
lines) are also shown. 

polynomial for the background. Centroid and width are 1.869 ±0.002 GeV/c 2 
and 0.012 ± 0.002 GeV/c 2 for D+ (Fig. 1(a)), and 1.970 ± 0.001 GeV/c 2 and 
0.010 ± 0.001 GeV/c 2 for D+ (Fig. 1(b)). The measured widths are in good 
agreement with the Monte Carlo predictions. The signal to noise ratio are 
1.0 ± 0.1 and 2.4 ± 0.4 for D + and D+, respectively. The reflected events are 
smoothly distributed across the invariant mass spectrum and are thus properly 
accounted for by the background polynomial fitting function. 

Residual structures induced by reflections in the Dalitz plot signal and side- 
band regions are a potential source of systematics of our amplitude analy- 
sis. The shape of the background in the Dalitz signal region is parametrized 
through a fit to the Dalitz plot of mass sidebands, which consists of a poly- 
nomial plus Breit-Wigner functions to account for any feed-through from res- 
onances. We verify that the reflections in the sidebands are well described by 
this simple background fit function and adequately represent the reflections 
in the signal regions. 



3 Branching ratio results 

The DCS D + — > K + tt + tt~ decay fraction is evaluated with respect to the 
CF D + -> K-n + n + mode, the SCS D+ -> K+it+it- with respect to the CF 
— ► K~K + n + mode. In Fig. 2 the mass distributions for the normalization 
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Figure 2. Invariant mass distributions of D + — ► ir + ir + and Z)+ — > K K + ir + 
normalization channels. 

channels are shown. The CF D + — > K~h + 'k + sample is selected with the 
same set of cuts as for the DCS channel, in order to minimize systematic 
effects; the signal yield is 32714 ± 184. The CF D+ -> K + K~>k + is obtained 
requiring for K~ a AWk > 2 similar to the AW n for the opposite-sign pion 
in the SCS selection; kinematical cuts are applied to remove the reflections in 
the K~K + Tr + mode from the CF D + -> K-n + ir + and from D* + -> D°ir + , 
followed by the D° decay to K~7r + . The — > K~K + it + signal yield is 
4033 ± 68. A fully coherent generation is used both for the channels under 
study and the normalization modes; the relative efficiencies are measured to 
be 



e(D + -> K+it+tt- 
e(D+ -»• K-n+n^ 

and 

t{D+ -> K + n + n- 



0.888 ±0.006 



1.106 ±0.009. (2) 



From the yields and efficiency ratios, we measure the following branching 
ratios, 



and 



0.0065 ±0.0008 ±0.0004 (3) 
= 0.127 ±0.007 ±0.014. (4) 



r (d+ -> ir+K-7r+) 

The first error reported is statistical and the second is systematic. 

To evaluate the systematic error we consider, as done in other FOCUS anal- 
yses [5,6], three contributions, which are added in quadrature to obtain the 
global systematic error: the split sample, fit variant and cut variant compo- 
nents. The split sample component takes into account the systematics intro- 
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duced by residual difference between data and Monte Carlo, due to a possible 
mismatch in the reproduction of the D momentum and the changing experi- 
mental conditions of the spectrometer during data collection. This component 
has been determined by splitting data into four independent subsamples, ac- 
cording to the D momentum range (high and low momentum) and the con- 
figuration of the vertex detector. The S-factor method from the Particle Data 
Group [7] was used to try to separate true systematic variations from statistical 
fluctuations. The branching ratio is evaluated for each of the four statistically 
independent subsamples and a scaled variance is calculated a (the errors are 
boosted when x 2 /{N — 1) > 1). The split sample variance is defined as the 
difference between the reported statistical variance and the scaled variance, if 
the scaled variance exceeds the statistical variance. Another possible source 
of systematics uncertainty is the fit variant. This component is computed by 
varying the fitting conditions on the whole data set. In our study fit variants 
include the background and signal shape, the bin size of the mass-distribution 
histogram and the Monte Carlo generation modeling. More precisely different 
degrees of the polynomial functions are used for the background parametriza- 
tion and two Gaussian peaks with the same mean but different widths are used 
for the signal, to take into account the different resolution in momentum of 
our spectrometer. The fully coherent generation for the decays under analysis 
is based on the results of the Dalitz plot study presented in this paper (see 
paragraph 4.3). To possible systematics effect in the branching ratio 

evaluation coming from our amplitude analysis we vary coefficients and phases 
returned by the Dalitz plot fit within their errors. The branching ratio values 
obtained by these variants are all a priori likely, therefore this uncertainty 
can be estimated by the standard deviation of the measurements. To inves- 
tigate the dependence of the branching ratio on the particular choice of cuts 
we calculate a cut variant error, analogously to the fit variant, by using the 
standard deviation of different branching ratio measurements obtained with 
several sets of cuts; this error is actually overestimated since the statistics of 
the cut samples are different. 

For the D + the dominant systematic effect comes from the cut variant, while 
for the the main source of uncertainty comes from the split sample compo- 
nent. Our results are consistent with previous determinations from the E687 [8] 
and E791 [9] experiments. Our measurements reduce the statistical errors by 
about a factor of 2 and 5 for the D + and D+, respectively. 



4 Dalitz plot analysis 

Events within ±2cr of the mass peak are used to perform the Dalitz plot 
analysis for D + and — > K + 7i + n^; the D + and Df Dalitz plots are shown 
in Fig. 3 a) and b), respectively. 
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Figure 3. a) D + and b) Dalitz plots in K + ir~ and 7r + 7r~ squared mass contri- 
butions. 

4-1 The decay amplitude 



For this analysis we use a formalization of the decay amplitude based on 
the simple isobar model, which is viewed as an effective description of the 
main dynamical features of the decay. More rigorous treatments, such as that 
based on the K-matrix formalism [10], are not viable for this analysis because 
of the limited statistics of the samples, along with the large number of free 
parameters necessary to account for the simultaneous presence of both n + n~ 
and K + 7t~ resonances. 

The decay matrix element for the analysis reported here is written as a co- 
herent sum of amplitudes corresponding to a constant term for the uniform 
direct three-body decay and to different resonant channels: 

M = a e i5 ° + Oj e ^ B(abc\r) (5) 
j 



where a, b and c label the final-state particles. B(abc\r) = B(a,b\r)S(a,c) 
where B(a, b\r) is the Breit-Wigner function 

B(a, b\r) = — (6) 

1 ' 1 ; M r 2 - M 2 ab -iYM r 1 ; 



and S(a, c) = 1 for a spin-0 resonance, S(a, c) = —2c ■ a for a spin-1 resonance 
and S(a, c) = 2(|c| |a|) 2 (3 cos 2 9* — 1) for a spin-2 state, a and c are the three- 



8 



momenta of particles a and c measured in the oh rest frame, and cos 6* = 
c ■ a/|c||a|. The momentum-dependent form factors Fp and F r represent the 
strong coupling at each decay vertex. For each resonance of mass M r and spin 
j we use a width 



where p is the decay three-momentum in the resonance rest frame and the 
subscript r denotes the on-shell values. The order of the particle labels is im- 
portant for defining the phase convention; here the first particle is the opposite- 
sign one. 

For the / (980) a simple single-channel Breit-Wigner is used. The mass and 
width values of 0.972 ± 0.002 GeV and 0.059 ± 0.004 GeV, respectively, are 
obtained from a fit to our FOCUS Dalitz plot of the Df — > 7r + 7r~7r + chan- 
nel, where the /o(980) is the dominant component. Given the questionable 
assumption of a single channel Breit-Wigner, where a coupled-channel one 
should be more properly used, we have checked that the final results do not 
change when a Flatte parametrization [11] is used. 

4-2 The likelihood function and fitting procedure 

Following our previous analyses [10,12,13] we perform a maximum likelihood 
fit to the Dalitz plots to measure the coefficients of the various decay am- 
plitudes as well as their relative phases. The probability density function is 
corrected for geometrical acceptance and reconstruction efficiency. The shape 
of the background in the signal region is parametrized through an incoherent 
sum of a polynomial function plus resonant Breit-Wigner components, which 
are used to fit the Dalitz plot of mass sidebands. The number of background 
events expected in the signal region is estimated through fits to the K + n + n^ 
mass spectrum. All background parameters are included as additional fit pa- 
rameters and tied to the results of the sideband fits through the inclusion in 
the likelihood of a x 2 penalty term derived from the covariance matrix of the 
sideband fit. The D + and samples are fitted with likelihood functions C 
consisting of signal and background probability densities. Checks for fitting 
procedure are made using Monte Carlo techniques and all biases are found to 
be small compared to the statistical errors. The systematic errors on our results 
are evaluated following the strategy already explained for the branching ratio 
measurements, defining a split-sample and a fit-variant component. The as- 
sumption that the shape of the background in the sideband is a good represen- 
tation of the background in the signal region potentially constitutes a source 




M r F?( P ) 

M ah F?(p r ) 



(7) 
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of systematic error. We take into account this effect in the fit-variant error, 
by varying the polynomial function degree and adding/removing the Breit- 
Wigner terms, which are introduced to take into account any feed-through 
from resonances in the background, and computing the standard deviation of 
the different results. 

In our Dalitz plot analysis we allow for the possibility of contributions from 
all known K + tt~ and 7r + 7r~ resonances [7] and from a flat non-resonant contri- 
bution accounting for the direct decay of the D mesons into three-body final 
states. The fit parameters are amplitude coefficients aj and phases 5j. 2 The 
general procedure, adopted for all the fits reported here, consists of several 
successive steps in order to eliminate contributions whose effects on our fit are 
marginal. Contributions are removed if their amplitude coefficients, ao and Oj 
of Eq. 5 are less than 3 a significant and the fit confidence level increases due 
to the decreased number of degrees of freedom in the fit. The fit confidence 
levels (C.L.) are evaluated with a x 2 estimator over a Dalitz plot with bin size 
adaptively chosen to maintain a minimum number of events in each bin. Once 
the minimal set of parameters is determined, addition of each single contribu- 
tion previously eliminated is reinstated to verify if the C.L. improves; in this 
case the contribution is added in the final set. 



4-3 Results 

The D + and Dalitz plots present a different event intensity distribution. 
The D + appears highly structured and dominated by p(770) and if* (892). The 
D+ channel, instead, presents a much more complex structure diffused over all 
the phase-space together with p(770) and if* (892), whose main characteristics 
are still visible. These general features are confirmed by our fits. 

In Table 1, fit fractions, 3 phases and coefficients of the various amplitudes 
describing the D + — > K + -n + -n~ decay are reported. The vector resonances 
p(770) and if* (892) account for about 90% of the D + decay fraction. Their 
relative phase difference, almost real, suggests a marginal role for final state 
interactions in this decay. The C.L. of our Dalitz plot fit is 9.2%. The three- 
lobe helicity structure visible in the D + Dalitz plot justifies the presence of 
the tensor if* (1430), at more than three sigma statistical significance, in the 
final set of resonances. The band of events at about (1 GeV/c 2 ) 2 in the m 2 +7r - 
mass combination indicates the presence of the scalar / (980) in the decay at a 

2 Coefficients and phases of the p(770), which is a dominant mode in both decays, 
are fixed to 1 and 0, respectively. 

3 The quoted fit fractions are defined as the ratio between the intensity for a single 
amplitude integrated over the Dalitz plot and that of the total amplitude with all 
the modes and interferences present. 
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X 2 = 13.6 
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Table 1 



Dalitz plot fit results for the D + — > if + 7r + 7r final state. 

four sigma significance level. The absence in the K + n~ combination of scalar 
1^q(1430), while the tensor ii"! (1430) is present, seems to us a bit suspicious; 
further investigations would required higher statistics. It is nevertheless inter- 
esting to observe that the DCS decay D + — > if + 7r + 7r~ is dominated by vector 
resonances with no major role of rescattering effects. The results of the fit on 
the two invariant mass squared projections m 2 Kn and are shown in Fig. 4. 

An amplitude analysis of this decay has been previously performed by the 
E791 experiment [9], which described the decay with two resonant channels, 
p(770)if + and if*(892)7r + , plus a uniform non-resonant component, each ac- 
counting for about 1/3 of the decay fraction; the non-resonant contribution 
seems to be better resolved in the present analysis in additional resonant chan- 
nels. We observe that the E791 p(770)/if *(892) relative phase difference 4 is 
about 0° while our measurement is close to 180°. We also perform a fit of 
our sample with the same components as E791. The results still confirm the 
disagreement. 

In Table 2, the fit fractions, phases and coefficients resulting from the 
fit are reported. As for the D + , the p(770) and if* (892) vector resonances 
represent the major contributions and cover about 60% of the decay fraction; 
their phase shift configuration is almost real for this decay as well. The de- 
scription of the event intensity all over the Dalitz plot requires three higher 
mass resonances, two in K + ix~ and one in 7r + 7r~, and a non- resonant term. 
More precisely the two K + n~ states are the scalar if,* (1430) and the vector 
if* (1410), which are the lowest mass resonances besides the if* (892), and the 
7r + 7r~ state is p(1450), which is the second vector state in the p series. They ac- 
count for an additional 30-40% resonant portion of the decay; a non-resonant 
contribution of about 15% completes the event description for this channel. 
The fit C.L. for the mixture of states selected is 5.5%. 

This solution satisfactory reproduces the main features of the decay, as indi- 
cated by the C.L., and shown in the two invariant mass squared projections 

4 This relative phase difference is free of phase convention ambiguity. 
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Table 2 

Dalitz plot fit results for the Df — » fr + 7r + 7r _ final state. 




m 2 K + n mass projection (GeV 2 /c 4 ) m 2 Jt + ^ mass projection (GeV 2 /c 4 ) 



Figure 4. The Dalitz plot projections and, superimposed, the fit results for D + . The 
background shape under the signal is also shown. 

m 2 Knr and m 2 n of the Dalitz plot in Fig. 5. However the absence of the /o(980) in 
the fit is a bit suspicious; an accumulation of events at (1 GeV/c 2 ) 2 tt + it~ mass 
squared, to some extent visually recognizable in the Dalitz plot, would indeed 
suggest its selection in the resonance final set. On the other hand we know 
that the isobar model is too naive to describe more complex decays dynamics, 
which intervenes in the presence of the K + n~ and 7r + 7r~ S-waves states. More 
rigorous treatments, such as that based on the K-matrix model [10], will be 
necessary at higher statistics. The results of the fit on the squared invariant 
masses m 2 KlT and are shown in Fig. 5. 

The dominant systematic uncertainties on coefficients and phases for both D + 
and decays come from the low/high momentum split and the fit variant 
component. 
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m 2 K + n - mass projection (GeV 2 /c 4 ) m 2 / n - mass projection (GeV 2 /c 4 ) 



Figure 5. The Dalitz plot projections and, superimposed, the fit results for Df. The 
background shape under the signal is also shown. 

5 Conclusions 



We have presented a study of the doubly and singly Cabibbo suppressed decays 
D + and — > K + tt + tt~. Our measurements of T(D + — » iC 1 ~7T + 7T~)/r(.D + — > 
ir-7r + 7r + ) = 0.0065 ± 0.0008 ± 0.0004 and T(D+ -> ^+7r + 7r-)/r( J D+ -> 
.K' + .R' _ 7r + ) = 0.127 ± 0.007 ± 0.014 improve the statistical accuracy by ap- 
proximately a factor of 2 and 5 with respect to previous determinations. In 
particular the comparison of (l/tan 4 # c ) x T(D£ cs )/T(D^ f ) = 2.60 ± 0.32 
with the FOCUS lifetime ratio of t(D+)/t(D°) = 2.538 ± 0.023 and the 
marginal role of FSI inferred by our Dalitz plot analysis of this decay sup- 
ports the interpretation that destructive interference between spectator am- 
plitudes with indistinguishable quarks in the CF D + final state is responsible 
for the lifetime difference between D + and D°. The FOCUS collaboration has 
already analysed the doubly Cabibbo suppressed decay of the neutral meson 
D° — > K + n~ [14]; the study of the DCS decay of the charged D + meson 
presented in this paper, free of any possible uncertainty due to D°D° mix- 
ing effects, provides complementary information. The amplitude analysis of 
D + and D£ —>■ K + tt + tt~ final states have also been performed, providing 
the first amplitudes and phases for the channel. We have discussed the 
major achievements and indicated possible improvements for a better decay 
dynamics interpretation at higher statistics. 
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